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摘要 : 采用 SSTKk-@9 模 型 对 冷却 条 件 下 超 临 界 压 力 CO, 在 水 平 管内 的 对 流 换 热 进行 了 数值 研究 , 分 析 了 流体 物性 、 热 流 密度 、 
直径 以 及 浮 升 力 等 对 其 在 拟 临界 点 附近 的 流动 换 热 特性 的 影响 ， 并 从 场 协同 的 角度 分 析 了 超 临界 压力 CO, 的 传 热机 理 。 结 果 表 
明 : 浮力 效应 使 流体 在 流动 截面 上 出 现 温 度 场 不 对 称 和 二 次 流 现象 ; 下 壁面 的 对 流 换 热 系数 比 上 壁面 先 达到 峰值 ， 但 换 热 系数 
小 于 上 壁面 ; 增 大 热流 密度 对 换 热 系数 的 影响 较 小 但 能 够 使 换 热 系数 的 峰值 向 入 口 段 移动 ; 增 大 热流 密度 和 增 大 直径 能 够 增强 
浮力 效应 对 流体 换 热 特性 的 影响 ， 场 协同 原理 可 以 解释 同一 截面 处 的 换 热 不 均匀 现象 。 
关键 词 ; 超 临界 压力 CO，; 对 流 换 热 ， 场 协同 原理 ; 浮力 效应 ;数值 模拟 
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A study on the cooling heat transfer mechanism for supercritical pressure 


C0, in horizontal tube 
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Abstract: The convective heat transfer of supercritical pressure CO, in horizontal tube under cooling condition is 
numerically investigated using SSTk-wturbulent model. The effects of thermal physical properties, heat flux, tube 
diameter and buoyancy on the heat transfer characteristics are discussed. The results show that the temperature is 
asymmetric and secondary flow is generated on the cross section due to the buoyancy effect. The peak value of heat 
transfer coefficient appears earlier on the bottom surface than on the top surface, but its value is smaller than that on 
the top surface. The increasing heat flux has little influence on the peak value of heat transfer coefficient, but affects 
its position seriously. The larger the heat flux and the tube diameter are, the more significant the buoyancy effect is. 
The non-uniformity of heat transfer on the same cross section can be explained by field synergy principle. 

Key words: supercritical pressure CO2>; convective heat transfer; field synergy principle; buoyancy effect; numerical 
simulation 
| 守 应 堆 、 太 阳 能 系统 、 制 冷 系统 等 领域 中 具有 广阔 的 
应 用 前 景 "*。 与 传统 常 物性 流体 相 比 ， 超 临界 CO， 
CO 流体 具有 无 毒 、 不 燃 、 性 能 稳定 以 及 价格 便 。 流体 在 临界 点 和 拟 临 界 点 附近 物性 变化 剧烈 ， 传 热 
宜 等 优点 ， 而 且 超 临界 CO, 流体 由 于 其 特殊 的 物性 ”规律 特殊 而 复杂 ， 因 此 超 临 界 CO, 流体 的 流动 换 热 
变化 ， 表 现 出 良好 的 流动 和 传 热 特 性 ， 因 此 在 核 反 特性 成 为 近年 来 的 研究 热点 ”。 
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I 研究 


Dang & Hihara 采用 实验 和 数值 模拟 方法 研 
究 了 直径 、 质 量 流 量 和 热流 密度 等 因素 对 超 临 界 压 
力 CO, 流体 在 水 平 管 中 的 对 流 换 热 特性 ， 并 且 假 设 
浮力 效应 忽略 不 计 。Dnu 等 人 "发 现 水 平 管 中 浮 力 效 
应 强化 了 超 临 界 压力 CO; 流体 在 拟 临界 点 附近 的 换 
热 强度 。 

过 增 元 等 


“提出 了 场 协 同 原理 ， 认 为 热源 强度 


不 仅 取决 于 流体 的 物性 用 


I 速 度 ， 还 取决 与 速度 和 尝 


1 


度 梯度 场 的 协同 性 。 当 速度 矢量 和 温 


度 梯度 矢量 的 


夹 角 小 于 900” 时 ， 夹 角 越 小 ， 


则 换 热 强度 越 好 。 


@ 7 

es- lo G) 
其 中 , 5 为 烙 性 引起 的 能 量 耗 散 , /4 为 基于 湾流 模型 
的 满 流 精度 。 


1.3 数值 方法 和 边 珊 条件 

采用 ANSYS CEFX ”软件 进行 数值 计算 ， 采 用 
NIST Standard Reference Database 23(REFPROP) 
Version 7 数据库 计算 超 临 界 CO, 的 热 物性 参数 ， 
并 编写 RGP 物性 文件 ， 导 入 CFX 软件 。 与 线性 插值 


本 文采 用 数值 方法 丰 


究 超 临界 压力 CO, 在 拟 临 


方法 相 比 ， 物 性 RGP 文件 能 够 更 精确 


地 反映 物性 的 


界 点 (ZT.=307.8 K，p=8 MPa) 附 近 的 对 流 换 热 特性 ， 


EE 点 分 析 热 


II 四 


\ 流 密度 、 直 径 和 浮 升力 对 局 部 对 流 换 热 


强度 的 影响 ,并 从 场 协 同 的 角度 分 析 超 临界 压力 CO。 


的 局 部 传 热机 到 


E， 以便 能 够 为 高 效 换 热 器 的 开发 和 


设计 提供 理论 基础 。 
1 数值 模拟 
1. 1 物理 模型 


由 于 超 临 界 压力 CO, 物性 变化 剧烈 ， 


为 了 准确 


7 计算 采用 压力 速度 耦合 算法 和 SSTK-w 模 

，SSTK-w 模 型 兼 具 有 大 -模型 的 健壮 性 和 独立 性 
ee 
制 方程 的 迭代 残 差 小 于 10* 时 ， 认 为 数值 计算 已 经 


条 件 ， 质量 流量 入 口 ， 压 力 出 口 ， 绝 热 入 
是 质量 流量 100-300 
kg/m2s， 压 力 8 MPa， 入 口 温度 340.15 人， 热流 密度 
35-45 kWym2， 入 口 雷 诺 数 约 为 3X10?。 


模拟 其 在 水 平 管内 的 局 部 对 流 换 热 特性 ， 建 立 三 维 


1.4 网 格 无 关 性 验证 和 数值 方法 验证 


模型 ， 


如 图 1 所 示 。 模 型 尺寸 : 直径 2、4、6 mm 的 


水 平 单 管 ， 
600 mm 冷却 段 组 成 ， 
近 充 分 发 展 ， 减 小 入 口 效 应 。 


长 度 840 mm， 由 240 mm 绝热 入 口 段 和 
以 便 在 冷却 段 开 始 处 使 流动 接 


网 格 由 ANSYS ICEM 生成 ， 径 向 网 格 加 密 ， 确 
保 近 壁面 y+<1， 以 满足 SST 湛 流 模型 对 y+ 的 要 求 ， 
网 格 质量 大 于 0.6， 满 足 要 求 。 网 格 无 关 性 验证 结果 
如 表 1 所 示 。 序 号 3 和 5 的 两 套 网 格 对 应 的 换 热 系 
数 的 相对 误差 为 0.2%， 考 虑 计算 时 间 和 求解 精度 ， 


Adiabatic 


选择 3 号 网 格 进行 模拟 。 为 了 验 训 


FE 数值 模拟 的 准确 


性 ， 将 模拟 结果 与 Dang & Hihara 


“的 实验 数据 进行 


-I 


= 六 | 全 
站 二 Ty 0=90° 
0=0° 


图 


1 物理 


模型 


Fig.1 Physical model 


1. 2 控制 方程 


直角 坐标 系 下 的 控 和 


| 方程 如 下 : 


对 比 。 模 型 尺寸 和 边界 条 件 均 与 文献 实验 工 况 一 致 。 
模拟 与 实验 结果 对 比如 图 2 所 示 。 模 拟 数 据 和 实验 
数据 的 最 大 相对 误差 为 13.5%, 说 明了 数值 方法 的 准 
确 性 和 可 靠 性 。 

表 1 网 格 无 关 性 验证 


Table 1 Mesh independence verification 


连续 性 方程 : 


Em pT (1) 


序号 。 网 格 数 y+ 1Wm2K7 ”误差 (%) 
1 156975 -0.743 10259.66 4.9 
2 268128 0.744 10188.78 4.2 
3 557403 0.594 9803.29 0.2 
4 697103 ”0.594 9803.29 0.2 
5 940615 0.595 9778.95 0 


10 上 5 了 
”Dang&Hihara 数据 


模拟 数据 ] 
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图 2 数值 方法 验证 六 


Fig. 2 Verification ofnumerical method 


2 场 协 同 原 理 


对 于 一 般 三 维 模型 ， 考 虑 稳 态 无 内 热源 情况 ， 
其 对 流 换 热 能 量 方程 的 积分 式 为 “: 


a 67 67 67 oT 67 
| PC| 2 一 -+Y 一 一 十 甸 kl 一 + 一 了 | 多 
0 Ox Oy Oz Ox Cz 


oT 
至 :二 大 = 一 | .三 
Oy dy 
其 中 ，D、c 和 大 分 别 为 流体 密度 、 
系数 ， 6 为 热 边 界 层 厚度 ，gqy 表示 壁面 热流 密 
忽略 轴 向 导热 ， 将 (4) 式 写成 矢量 形式 : 


(4) 


pe (U :VT)dy= gq, (5) 

引入 无 因 次 变量 : 

二 去 VT 一 Jy 

U=—, V7T=—, = 一 (6) 
i [= 


其 中 ，wb 和 分 别 为 流体 平均 速度 和 平均 温度 ，T 

为 壁面 温度 。 

将 (6) 式 代入 (5) 式 并 整理 可 得 无 因 次 关系 式 : 
mW (7) 

(T, -下 ) 有 


其 中 , h 为 对 流 换 热 系数 ,雷诺 数 Re = xd& ， 普 朗 特 
Hu 


jiRe, Pn, (UvT)dy= 


数 Py = 全 ,下 标 b 表示 主流 体 ,下 标 w 表示 壁面 ， 
上 .V7 可 以 表示 为 : 
U-.v7=|U|-IvTlcosp (8) 


其 中 ，p 为 速度 矢量 与 温度 梯度 矢量 之 间 的 来 角 ， 即 
场 协同 角 。 


广 im”Avi (大 甘 日 工 
CNINAAIVD TFANDT! 


式 (7) 表 明 在 变 物性 情况 下 ，Nw, 与 流体 的 流动 、 
物性 和 通道 直径 有 关 ， 还 与 速度 场 和 温度 梯度 场 的 


协同 性 有 关 , 影响 因素 众多 ,比较 复杂 。 因 此 将 式 (7) 

改写 成 : 

fp 本 (9) 
计算 结果 与 讨论 


3. 1 热流 密度 对 流动 换 热 的 影响 

为 了 研究 热流 密度 对 超 临 界 压 力 CO, 流动 换 热 
特性 的 影响 ， 采 用 直径 为 6 mm 热流 密度 为 35、45 
kW， 压力 8 MPa， 入 口 温 度 340.15 K。 

图 3(a) 给 出 了 不 同 热 流 密 度 条 件 下 , 流体 平均 温 
度 Th 和 局 部 壁面 温度 7, 沿 程 分 布 。 图 3(a) 表 明 在 定 
热流 条 件 下 ， 流 体温 度 和 壁面 温度 沿 流向 逐渐 降低 
到 拟 临 界 温度 以 下 ; 上 壁面 温度 比 下 壁面 温度 下 降 
快 ， 这 是 由 于 超 临 界 压 力 CO; 物性 变化 剧烈 ， 在 水 


A te ee 使 得 上 下 
壁 温 分 布 并 不 对 称 ; 当 热流 密度 增 大 时 ， 流 体 的 冷 


ne tt et 
说 明 增 大 热流 密度 能 够 增强 浮 升力 对 超 临 界 压力 
CO, 流动 换 热 的 影响 。 

图 3(b) 给 出 了 对 流 换 热 系数 有 沿 程 分 布 ， 换 热 
系数 调 流 向 先 增 大 到 峰值 后 减 小 ， 这 与 超 临 界 CO， 
的 定 压 比 热 随 温 度 的 变化 趋势 一 致 。 下 壁面 换 热 系 
数 比 上 壁面 先 达 到 峰值 ， 但 其 数值 明显 小 于 上 壁面 ， 
在 冷却 条 件 下 ， 靠 近 壁 面 的 流体 先 被 冷却 ， 由 于 浮 
力 影 响 ， 温 度 高 而 密度 小 的 流体 向 上 移动 ， 导 和 臻 上 
壁面 温度 与 主流 体 之 间 的 温差 较 小 ， 因 此 上 壁面 具 
有 更 高 的 换 热 系数 。 当 热流 密度 增 大 时 ， 局 部 换 热 
系数 的 峰值 几乎 没有 变化 ， 但 达到 峰值 的 位 置 提前 ， 
这 是 由 于 增 大 热流 密度 使 流体 温度 更 快 降低 到 拟 临 
界 温 度 ， 在 拟 临 界 温 度 附 近 ， 换 热 系 数 达 到 峰值 。 
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图 3 不 同 热 流下 (a) 温 度 和 (b) 对 流 换 热 系数 沿 程 分 布 
Fig. 3 Distributions of (a) temperature and (b) heat transfer 


coefficient along the tube at different heat fluxes 

图 4 给 出 了 不 同 热流 密度 下 局 部 Nem 分 布 情况 ， 
随 着 热流 密度 的 增 大 ，Nu 的 峰值 略 有 提高 ， 说 明 增 
大 热流 密度 有 利于 增强 对 流 换 热 强度 。 图 5 给 出 了 
该 工 况 下 局 部 场 协同 角 的 分 布 ， 场 协同 角 沿 程 分 布 
不 均匀 ， 近 上 壁面 处 的 场 协同 角 小 于 下 壁面 ， 说 明 
近 上 壁面 处 速度 和 温度 梯度 的 协同 性 比 下 壁面 好 ， 


gq.=-45 kW/m” 
top—®—bottom 
G39 kW/m” 
一 全 top 一 Y 一 bottom 


区 | 


5 不 同 热 流 密度 下 场 协 同 角 沿 程 分 布 
Fig. $ Distribution of f along the tube at different heat fluxes 
3. 2 浮 升力 对 流动 换 热 的 影响 

为 了 更 加 直观 的 了 解 浮 升力 对 超 临 界 压力 CO， 
流动 换 热 特性 的 影响 , 在 不 考虑 重力 (ge=0) 情 况 下 进 
行 了 模拟 ， 并 与 考虑 重力 (gj=-9.81 m/s ) 工 况 进行 对 
比 。 

图 6 给 出 了 重力 对 Nuo 沿 程 分 布 的 影响 。 从 图 中 
可 以 看 出 : 考虑 重力 时 的 局 部 对 流 换 热 强度 略 高 于 


因此 具有 更 好 的 换 热 特性 ;在 其 它 条 件 相 同时 ， 同 
一 截面 处 ， 上 下 近 壁 面 处 的 场 协同 角 之 间 的 差距 随 
着 热流 密度 的 增 大 而 增 大 ， 说 明 增 大 热流 密度 使 浮 
力 效 应 引起 的 对 流 换 热 场 的 不 均匀 性 更 加 明显 。 
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图 4 不 同 热流 密度 下 Nuo 沿 程 分 布 
Fig. 4 Distribution of Nu;, along the tube at different heat fluxes 


不 考虑 重力 时 的 工 况 ， 说 明 重 力 引 起 的 浮力 效应 能 
够 增强 换 热 强度 ， 但 主要 表现 在 冷却 的 后 半 段 ， 即 
主流 体温 度 小 于 拟 临界 温度 的 区 域 ， 随 着 热流 密度 
的 增 大 ， 重 力 引 起 的 浮力 效应 更 加 明显 。 图 7 给 出 
了 对 应 工 况 下 的 局 部 场 协 同 角 ， 可 以 看 出 ， 不 考虑 
E 力 时 的 场 协 同 角 小 于 考虑 重力 时 的 工 况 ， 因 此 对 
于 物性 剧烈 变化 的 超 临界 流体 ， 场 协同 原理 的 应 用 
不 能 简单 归 因 为 场 协同 角 的 减 小 ， 应 综合 考虑 多 种 
影响 因素 。 
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图 6 有 无 重力 条 件 下 Nam 沿 程 分 布 
Fig. 6 Distribution of Nzu along the tube with or without gravity 
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图 7 有 无 重力 条 件 下 场 协 同 角 沿 程 分 布 
Fig. 7 Distribution of f along the tube with and without gravity 


图 8 给 出 了 有 无 重力 条 件 下 ， 当 量 热源 的 数值 
jpe,(U .v7 ) 的 沿 程 分 布 。 考 虑 重力 时 的 当量 热源 明 


显 大 于 不 考虑 重力 时 的 工 况 ， 而 且 考 虑 重力 时 当量 
热源 沿 程 分 布 不 均匀 ， 当 量 热源 随 着 热流 密度 的 增 
大 而 增 大 。 根 据 式 (5) 和 (7)， 增 大 当量 热源 能 够 增强 
换 热 强 度 ， 当 量 热 源 的 数值 不 仅 与 速度 、 温 度 以 及 
两 者 的 夹 角 有 关 ， 还 与 流体 的 物性 有 关 ， 由 于 超 临 
界 CO, 流体 的 物性 随 温度 变化 剧烈 ， 因 此 物性 的 剧 
烈 变 化 成 为 其 流动 换 热 特 性 的 主要 影响 因素 。 
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图 8 有 无 重力 条 件 下 -oc, (UV .V7) 沿 程 分 布 
Fig. 8 Distribution of -pe, (U :VT ) along the tube with and 
without gravity 

3. 3 直径 对 流动 换 热 的 影响 

为 了 研究 直径 对 超 临 界 压力 CO; 流动 换 热 特性 
的 影响 ， 采 用 直径 为 2、4、6 mm 热流 密度 为 45 
kW/m ， 压 力 8 MPa， 入 口 温度 340.15K， 入 口 雷诺 
数 均 为 3X10 。 


图 9 给 出 了 不 同 直径 对 应 的 Nu 沿 程 分 布 ， 在 
Na 达到 峰值 之 前 ，Nuo 随 着 直径 的 增加 而 增加 ， 达 
到 峰值 之 后 ，Nu 逐渐 减 小 ， 此 时 直径 大 小 对 其 影响 
较 小 。 该 工 况 对 应 的 场 协 同 角 分 布 如 图 10 所 示 ， 由 
图 10 可 见 ， 流 体温 度 达 到 拟 临 界 温度 之 前 ， 场 协同 
角 随 着 直径 的 增 大 而 明显 减 小 ， 即 在 其 他 条 件 相同 
时 ， 增 大 管 径 能 够 改善 速度 场 和 温度 梯度 场 的 协同 
性 ， 增 强 换 热 强度 。 当 流体 温度 小 于 拟 临 界 温度 时 ， 
管 径 对 场 协同 角 的 影响 较 小 ， 同 时 对 传 热 系 数 的 影 
响 也 比较 小 。 
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图 9 不 同 直径 对 应 的 Nm 沿 程 分 布 
Fig. 9 Distribution of Nu, along the tube with different diameters 
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图 10 不 同 直径 对 应 的 场 协 同 角 沿 程 分 布 


Fig. 10 Distribution of f along the tube with different diameters 

图 11 给 出 了 超 临界 CO; 在 直径 2 mm 和 6 mm 
的 水 平 管内 冷却 时 ， 一 390 mm 截面 上 流体 温度 、 径 
向 速度 和 场 协同 角 的 分 布 。 截 面 上 流体 温度 出 现 不 
对 称 分 布 ， 温 度 高 而 密度 小 的 流体 向 上 聚集 ， 温 度 
低 而 密度 大 的 流体 主要 集中 在 底部 ;重力 引起 的 译 
力 效 应 使 得 截面 上 形成 二 次 流 ， 增 加 了 流体 的 混合 
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程度 ;， 截面 上 场 协 同 角 也 呈 不 均匀 分 布 ， 上 部 的 场 
协同 角 小 于 底部 的 场 协 同 角 ; 随 着 直径 的 减 小 ， 截 
面 上 温度 分 布 和 场 协 同 角 分 布 更 加 均匀 ， 二 次 流 强 
度 减 小 。 
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Fig. 11 Distributions oftemperature, velocity and field synergy 


angle on cross section 
4 结论 


本 文 利用 数值 方法 研究 了 冷却 条 件 下 超 临 界 压 
力 CO;, 在 水 平 管内 的 对 流 换 热 特性 ,主要 结论 如 下 : 

(1) 浮 力 效应 使 截面 上 流体 温度 场 分 布 不 对 称 ， 
同一 截面 处 上 壁面 温度 大 于 下 壁面 ， 近 壁面 流体 先 
被 冷却 ， 产 生 较 大 密度 差 ， 形 成 二 次 流 。 

(2) 下 壁面 的 对 流 换 热 系数 比 上 壁面 先 达 到 峰 
值 ， 但 换 热 系数 小 于 上 壁面 ， 截 面 顶部 的 场 协 同和 角 
小 于 底部 。 

(3) 热 流 密 度 的 增 大 使 换 热 系 数 的 峰值 向 入 口 段 
移动 ， 增 大 热流 密度 和 直径 能 够 增强 浮力 效应 对 流 
体 换 热 特性 的 影响 。 

(4) 与 g=0 相 比 ， 考 虑 重力 (@=-9.81 m/s”) 时 虽然 
场 协同 角 较 大 ， 但 当量 热源 也 比较 大 ， 因 此 换 热 强 
度 大 于 ge=0 时 的 工 况 。 
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